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ABSTRACT 

We present the most sensitive ultraviolet observations of Supernova 1987A to date. Imaging spec- 
troscopy from the Hubble Space Teles cope-Gosvuic Origins Spectrograph shows many narrow (Au ~ 300 
km s~ 4 ) emission lines from the circumstellar ring, broad (Aw ~ 10 - 20 x 10 3 km s _1 ) emission lines 
from the reverse shock, and ultraviolet continuum emission. The high signal-to-noise (> 40 per reso- 
lution element) broad Lya emission is excited by soft X-ray and EUV heating of mostly neutral gas in 
the circumstellar ring and outer supernova debris. The ultraviolet continuum at A > 1350 A can be 
explained by H I 2-photon (2s 2 Si/ 2 - Is 2 Si/ 2 ) emission from the same region. We confirm our earlier, 
tentative detection of N V A 1240 emission from the reverse shock and present the first detections of 
broad He II A1640, C IV A1550, and N IV] A1486 emission lines from the reverse shock. The helium 
abundance in the high-velocity material is He/H = 0.14 ± 0.06. The N V/Ha line ratio requires 
partial ion-electron equilibration (T e /T p w 0.14 - 0.35). We find that the N/C abundance ratio in the 
gas crossing the reverse shock is significantly higher than that in the circumstellar ring, a result that 
may be attributed to chemical stratification in the outer envelope of the supernova progenitor. The 
N/C abundance ratio may have been stratified prior to the ring expulsion, or this result may indicate 
continued CNO processing in the progenitor subsequent to the expulsion of the circumstellar ring. 
Subject headings: supernovae: individual (SN 1987A) — shock waves — circumstellar matter 



1. INTRODUCTION 

Borkowski, Blondin, & McCray (1997) predicted that 
the spectrum of SN 1987A should display very broad 
(Av ~ ± 12,000 km s~ 4 ) emission lines of Lya, Ha, 
N V A1240, and He II A1640, produced where the 
freely expanding supernova debris crosses a reverse 
shock located inside the equatorial circumstellar ring. 
In September 1997, using the Hubble Space Telescope- 
Space Telescope Imag ing Spectrograph (HST-STIS), 
ISonneborn et al.l (I1998D d etected broad Lya emission, 
and iMichael et al.l (|1998h showed how observations of 
this emission can be used to map the shape of the re- 
verse shock and the flux of H I atom s crossing it. Michael 
et al (2003) and lHeng et alj (|2006l ) analyzed subsequent 
(February 1999 - October 2002) STIS observations of 
both Lya and Ha to map the increasing flux of H I atoms 
across the reverse shock. 

The evolution of the broad Ha emission can also be 
tracked with ground-based telescopes, notably in Febru- 
ary 2005 with the Magellan telescope (| Smith et al.ll2005h 
and from December 2000 to January 2009 with the Very 
Large Telescope (Fransson et al 2011). Smith et al. 
pointed out that the Ha emission from the reverse shock 
could be suppressed due to preionization of hydrogen in 
the supernova debris by soft X-rays and EUV radiation 
from the rapidly brightening shock interaction with the 
inner circumstellar ring. Extrapolating the X-ray light 

Based on observations made with the NASA/ESA 
Hubble Space Telescope, obtained from the data archive at the 
Space Telescope Science Institute. STScI is operated by the As- 
sociation of Universities for Research in Astronomy, Inc. under 
NASA contract NAS 5-26555. 



curve, they predicted that this preionization would cause 
the broad Ha emission to vanish by 2012. Ho wever, since 
2005 , the X-ray light curve has leveled off (|Park et al.l 
2011), so the preionization effect should be less than that 
p redicted by Smith e t al. 

iFrance et al.l (|2010l ) analyzed the most recent (31 Jan- 
uary 2010) STIS observations of SN 1987A. They pointed 
out that the observed ratio of Lya/Ha photon fluxes 
was > 30 at large, negative velocities, much greater 
than the value « 5 that would be expected for hydrogen 
atoms excited as they cross the reverse shock. Moreover, 
the spatially resolved Lya line profile differed dramat- 
ically from that of Ha. France et al. proposed that 
the broad Lya emission is dominated by Lya emission 
from the nearly stationary equatorial ring that has been 
resonantly scattered by hydrogen atoms in the expand- 
ing supernova debris. France et al. also noticed a faint 
glow at wavelengths ranging from about 1260 - 1290 A, 
which they attributed to the N V A 1240 emission pre- 
dicted by Borkowski et al. They suggested that a critical 
test of that hypothesis would be detection of broad emis- 
sion from C IVA1550. 

Here we describe far-ultraviolet (1140 - 1780 A) spec- 
tra of SN 1987A obtained with the Cosmic Origins Spec- 
trograph (COS) that was installed on HST in May 2009. 
With COS, we have measured the profiles of Lya and 
N V A1240 with signal-to-noise ratios (S/N) far superior 
to the previous STIS spectra. We have also detected the 
broad He II A1640 emission predicted bv lBorkowski et all 
(jl997t ) and the broad C IV A1550 emission predicted by 
France et al (2010). We find that while the He/H abun- 
dance ratio in the outer ejecta is consistent with that 
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Fig. 1. — (top) H ST- ACS BVR image with a representative over- 
lay of the 2.5" diameter COS aperture (cyan), (bottom) HST-COS 
near-UV pointing verification image for the G130M (2011 February 
11) visit. Due to bright limits, this observation was made using the 
MIRRORB imaging configuration, which creates the image dou- 
bling. Because COS is a slitlcss spectrograph, light from outside 
the nominal 2.5" can enter the system. Star 3, separated from the 
center of the circumstellar ring by fa 1.7" is seen at the bottom. In 
order to match the spectral image presented in Figure 2, we display 
the coordinates of this image with east to the right of north, op- 
posite from the standard convention. The spectroscopic dispersion 
axis is labeled "A" and the cross-dispersion axis is labeled "y" . 

of the circumstellar ring, the N/C ratio is enhanced by 
more than a factor of two over the ring abundances. 

2. HST-COS OBSERVATIONS, IMAGING SPECTROSCOPY, 
AND DATA REDUCTION 



SN 1987A was observed with the medium resolution 
far-UV modes of HST-COS (G130M and G160M) on 
2011 February 11 and March 14 for a total of 7 spacecraft 
orbits (18555 s; Table 1) as part of the Supernova 1987A 
INTensive Study (SAINTS - G012241; PI - R. Kirshner). 
A description of the COS instrument a nd on-orbit per- 
forma nce characteristics can be found infOstcrma n et all 
(|2011l ) . All observations were approximately centered on 
the SN 1987A circumstellar ring (R.A. = 05 h 35 m 28.07 s , 
Dec. = -69°16'10.8"; J2000) and COS performed an off- 
set imaging target acquisition from a reference star w 8" 
away. A NUV image was obtained with the MIRRORB 
imaging mode after the slew to the science pointing in 
order to verify that the circumstellar ring was in the pri- 
mary science aperture (PSA). Figure 1 (top) shows the 
COS aperture location on an optical image of SN 1987A, 
The MIRRORB configuration introduces optical distor- 
tions into the image, but a first-order analysis indicates 
that the entire circumstellar ring was inside the 2.5" di- 
ameter PSA (Figure 1, bottom). 

The COS position angle was chosen to maximize the 
spatial separation on the COS detector between Star 3 
and the supernova emission we aim to study. The G130M 
observations were made at a position angle of s» —75°, 
while the G160M observations were acquired with a po- 
sition angle of « —45°. Light from objects outside the 
nominal 1.25" COS aperture radius can be recorded with 
the science spectru m in crowded fields. Star 3 is a far- 
UV-b right Be star (jGilmozzi et al.lll987t iWalborn et~aTl 
Il993f l separated from the center of the SN 1987A ring by 
~ 1.7". Stellar contamination could compromise super- 
nova data quality if the two objects overlapped on the 
spectroscopic detector. Thus, we used the imaging spec- 
troscopic capability of COS to keep SN 1987A centered 
while placing Star 3 at the bottom of the microchanncl 
plate (MCP) detector. 

We complemented this observing strategy with post- 
processing techniques and custom spectral extractions to 
maximize the spatial resolution of the instrument in the 
cross-dispersion direction. Due to long-term exposure 
of the COS MCP to geocoronal Lya and hot-star spec- 
tra, the number of electrons generated by each incident 
photon at a given detector location has been decreasing. 
This is manifest as lower pulse heights observed by the 
detector electronicsQ. The exact form of the pulse-height 
evolution is a complicated function of the dispersion and 
cross-dispersion arrival location of the incident photon, 
but to first order we can approximate this as a pulse- 
height dependent mislocation of the photon y-location. 
The approximate form of this mislocation is -0.47 y pix- 
els per pulse height bin drop. We have reprocessed the 
COS observations of SN 1987A, including a correction 
for the cross-dispersion misplacement, producing more 
accurate photon arrival positions. This additional data 
reduction step improves the far-UV angular resolution of 
the instrument by 20 - 30% across the bandpass. 

Star 3 can be seen at the bottom of the two- 
dimensional spectrogram shown in Figure 2. Star 3 is 
well isolated at the bottom of the detector. The angu- 
lar resolution in the reprocessed two-dimensional data 

2 We refer the reader to the cycle 19 

COS Instrument Handbook for more details: 

http : / /www. stsci . edu/hst/cos/documents/handbooks/current/co 
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is « 0.8". This imaging capability allows us to make 
custom spectral extractions at three spatially resolved 
locations: the "WNW" region at the top of the detec- 
tor (approximately centered on the hotspots at: R.A. 
= 05 h 35 m 27.98 s , Dec. = -69<T6'10.5"; J2000), the 
"ESE" region towards Star 3 (approximately the tele- 
scope pointing coordinates), and a separate extraction 
of Star 3 itself. Star 3 contributes negligible flux to the 
one-dimensional WNW spectrum and less than 15% of 
the flux in the ESE extraction. The spatially resolved 
spectra will be analyzed in a future work, and we focus 
on the combined WNW + ESE spectra here. 

The custom data extractions were then reprocessed 
with the COS calibration pipeline, CALCOS v2.13.6, and 
combined with the custom IDL c oadd ition procedure de - 
scribed bv lDanforth et al.l (pOlOl ) and lShull etHI (pOlOt ). 
The coaddition routine interpolates all detector segments 
and grating settings onto a common wavelength grid, 
and makes a correction for the detector QE-enhancement 
grid. No correction for the detector hex pattern is per- 
formed. Data were obtained in four central wavelength 
settings in each far-UV grating mode (A1291, 1300, 1309, 
and 1318 with G130M and A1577, 1589, 1600, and 1611 
with G160M) at the default focal-plane split position. 
The instrumental configurations are summarized in Ta- 
ble 1. Observations at multiple wavelength settings pro- 
vide continuous spectral coverage over the 1136 - 1782 A 
bandpass and minimize the residual fixed pattern noise 
from the detector grid wires and the MCP pores. The 
point source resolving power of the medium resolution 
COS far-UV modes is R = AA/A « 18,000 (Av = 17 
km s _1 ); however, the filled-aperture resolving power is 
R ~ 1500 (Av ~ 200 km s" 1 ; France et al. 2009). Multi- 
ple point sources within the aperture (e.g., the multiple 
hotspots along the SN 1987A circumstellar ring) produce 
spectral resolution that is similar to the extended source 
response. The point source flux calibration of COS is 
accurate to better than w 3%; however, the extended 
source and custom processing limit the absolute flux cal- 
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Fig. 2. — Two-dimensional spectrogram of the G130M segment B 
data. Reverse shock emission is labeled in yellow, hotspot emission 
is labeled in blue, and interstellar absorption features are marked 
in green. 
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Fig. 3. — Full far-UV spectrum of the combined WNW and ESE 
regions. Broad emission from the SN 1987A reverse shock, narrow 
emission lines attributable to circumstellar hotspots, as well as 
underlying continuum emission are observed. Representative 1-a 
error bars (a combination of photon statistics and flux calibration 
uncertainties) are shown in red. 



ibration of the SN 1987A far-UV spectroscopic data to 
~ 10%. 

3. DATA ANALYSIS 
3.1. Emission Spectrum and Interstellar Corrections 

The extracted emission spectrum of SN 1987A is shown 
in Figure 3. The spectra presented here are shown in 
photon units as these are the natural units for line flux 
comparisons. There is a wealth of narrow and broad 
emission and absorption features in the spectrum as 
well as a strong underlying continuum. The "narrow" 
(FWHM » 200 - 300 km s" 1 ) resolved emission fea- 
tures are attributed to the circumstellar ring. The cir- 
cumstellar ring spectrum is dominated by the hotspots 
which originate at shock interfaces where the forward 
blast wave encounters regions of high density in the cir- 
cumstellar ring (Lawrence et al. 2000). A discussion 
of the ring-ejecta interaction is given by Larsson et al. 
(2011). Most of the circumstellar ring emission lines have 
been observed in previous iJST-STIS (e.g., Pun et al. 
2002) and IUE (Lundqvist & Fransson 1996) observa- 
tions, including N V A1239, 1243, C II A1334, 1335, Si IV 
A1394, 1403, the O IV] A1400 multiplet, N IV] A1483, 
1486, C IV A1548, 1550, [Ne IV] A1602, He II A1640, 
O III] A1661, 1666, and the N III] A1750 multiplet. Our 
deep COS observations enable first detections of O V 
A1371 and Si II A1526, 1533, although the ground state 
transition of Si II is blocked by the interstellar medium. 
We also tentatively detect the coronal [Si VIII] doublet 
A1441, 1446. These features will be discussed in detail 
in a future work, though we will make use of the circum- 
stellar ring line fluxes where relevant. 

The broad emission features observed in the spec- 
trum of SN 1987A arise as atoms and ions cross the 
reverse shock front ([Michael et al.l 119981 ) , located just 
interior to the circumstellar ring. H I Lya (and the 
corresponding Ha ) emission from the reverse shock has 
been well studied ((Michael et alj|2003t iSmith et al.ll2005t 
iHeng et aT1l2006t iFrance et al.ll2010l) 7and we display the 
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TABLE 1 
SN 1987A COS observing log. 



Central Wavelengths FP-POS T exp (s) 

1291,1300,1309,1218 3 7884 

1577,1589,1600,1611 3 10671 



Date COS Mode 



2011 February 11 G130M 
2011 March 14 G160M 
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Fig. 4. — RS Lya and hotspot N V in the WNW extraction. 
Hotspot emission lines are marked in blue and interstellar absorp- 
tion components are labeled in green. The purple line shows a fit to 
the Lya line core for logio(N(H)) = 21.43 (± 0.02) cm -2 . Dashed 
orange lines are plotted as velocity references. Representative 1-a 
error bars (a combination of photon statistics and flux calibration 
uncertainties) are shown in red. 

broad (Av ~ —13,000 — 1-7,000 km s -1 ) Lya emission 
seen in the COS observations in Figures 4 and 5. The 
second broad feature to the red of Lya in Figure 3 is 
consistent wi th the putati v e N V reverse shock emission 
described bv lFrance et al.l (|2010l ) and has a velocity dis- 
tribution and total flux that are approximately consistent 
with those predicted by Borkowski et al. (1997); see also 
§4.3.1. In light of these results and the large relative off- 
set from the observed Lya profile, we claim that reverse 
shock N V emission is unambiguously detected in our 
data (Figure 5). While at earlier times, low-ionization 
UV line emission from SN 1987A has been a ttributed to 
the debris in the core ( Jcrks trand et al.ll2011l) , the core is 
enshrouded in dust (jMatsuura et al.ll2011l ) and is mostly 
opaque to UV radiation. Therefore, the reverse shock is 
the most likely formation site for high-velocity UV line 
photons observed towards SN 1987A. Following similar 
line-identification arguments as those for N V, we also 
detect reverse shock emission from N IV] A1486, C IV 
A1550, and He II A1640 for the first time. These features 
are marked in Figure 6. 

Figure 4 also shows an example of the broad and 
narrow interstellar absorption features imposed on the 
spectrum of SN 1987A. This figure shows Galactic and 
Magellanic Cloud absorption from Si II A1190, 1193, 
the N I A1200 multiplet, Si III A1206, and H I Lya. 
Other interstellar absorbers include S II A1250, 1253, 
Si II A1260, C II A1334, Si II A1526, Fe II A1608, 
and Al II A1670. We fit the interstellar neutral hydro- 
gen absorption in the WNW spectrum as this provides 
the best measure of the foreground absorption, finding 



logi (N(H)) = 21.43 ± 0.02 cm' 2 . For subsequent anal- 
ysis of the Lya and N V emission, we divide the data by 
a model of this emission (shown in purple in Figure 4). 
This procedure renders unusable the inner ± 6 A of the 
Lya line profile, which we remove from the data. 

There are several estimates of the interstellar redden- 
ing tow ards the region of the LMC in which SN 1987A 
resides dWalker fe Sunt zcfi 1990; iFitzpatrick fe Walbornl 
H990HScuderi et alJll996tlMichael et al.ll2003fl. We adopt 
a Mil ky Way (Ry = 3.1) extinction curve (jCardelli etaLl 
1989) with E{B — V) = 0.19 as this curve approximates 
the average extinction correction favored by the stud- 
ies cited above. This curve also closely approximates 
the far-U V properties o f the average LMC curve pre- 
sented by iGordon et ail (120031) for Ay = 0.6. Assum- 
ing E(B - V) = 0.17 (jMichael et"aT1[200l . would lower 
the fluxes de r ived i n subsequent sections by ~ 15%. 
iScuderi et al.1 (|1996h present a comprehensive study of 
the extinction towards Star 2, and the use of this curve 
would result in fluxes ~ 20 - 30% larger than those de- 
rived below. None of the possible extinction curves sig- 
nificantly alters the shape of the continuum (§4.1.1) or 
the far-UV line ratios (§4.3.3). 

The bright continuum (F\ (observed) ~ 1 x 10~ 4 pho- 
tons cm -2 s _1 A -1 ) underlying the hotspot and reverse 
shock emission was suggested in coadded STIS G140L 
spectra from 2010, cospatial with the circumstellar ring. 
However, the significance of the detection was low. We 
propose that this continuum is mostly H I 2-photon emis- 
sion (§4.1.1). 

3.2. Binned Spectrum for Broad Line & Continuum 
Analysis 

Narrow emission lines, as well as the narrow interstellar 
absorption components, complicate measurements of the 
continuum and broad emission from the reverse shock. In 
order to remove confusion from these narrow lines during 
the analysis of the broad spectral features, we created 
binned spectra of the WNW and ENE extractions. These 
binned spectra were centered on hand-chosen regions of 
the data at 1 - 5 A intervals and were free of narrow 
features. The flux and error of these binned data points 
were taken to be the average flux and standard deviation 
of a 0.8 A region centered on the selected wavelength. 
These binned data are shown as diamonds overplotted 
on the spectra in Figures 5-8. 

We quantify the total reverse shock flux by integrating 
the binned spectra over wavelength intervals correspond- 
ing to the velocity ranges of interest for a given feature. 
The velocity ranges were chosen to cover the maximum 
extent of a given line without significantly overlapping 
with other reverse shock emission features. The inner 
approximately 3000 km s _1 were not included to avoid 
strong hotspot emission lines located at ~ km s _1 . 
In most cases the hot spot lines are doublets, therefore a 
velocity interval larger than the nominal 300 km s" 1 nar- 
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Fig. 5. — Velocity profile of the ISM corrected (dust and neutral 
hydrogen), continuum subtracted RS Lya and N V in the com- 
bined WNW + ESE spectrum. The binned spectrum described 
in §3.2 are shown overplotted in gray. The Lyo emission extends 

from —13000 (-7000 km s — 1 and emission to the red of Lya is 

attributed to N V. 

row line-width was excluded. We were unable to measure 
two reverse shock components due to velocity blending. 
The blue component of the N V line is lost under the 
Lya line, and the weak red wing of the N IV] emission 
is overwhelmed by the stronger blue C IV emission. The 
velocity and wavelength intervals are given in Table 2 
with the summed line fluxes. 

4. RESULTS 

4.1. Hydrogen: 2-photon Emission and Broad Lya 
4.1.1. H I Continuum 

Figure 3 shows that the emission and absorption lines 
are superposed upon a far-UV continuum with F\ > 10~ 4 
photons cm~ 2 s _1 A -1 . This continuum was marginally 
detected in previous STIS G140L spectra, with the emis- 
sion being concentrated in the ring plane. It is detected 
in our COS observations at high significance, (Flux / 
Error) > 8 per spectral resolution element from 1360 - 
1540 A. 

We can rule out the possibility that this emission is 
uncorrected scattered light from Star 3. As described 
in Section 2, care was taken to isolate the emission from 
Star 3 at the bottom of the MCP detector (Figure 2). At 
the observed continuum levels, the scattered light from 
Star 3 would have to be over twice as bright outside of 
the stellar extraction region as within it, and relatively 
constant with angular separation from the stellar spec- 
trum. Additionally, the observed continuum of Star 3 
has a different spectral shape than the supernova con- 
tinuum and there is no evidence for stellar features in 
the supernova spectrum (such as photospheric and wind 
absorption lines of C III A1175 and Si IV A1394, 1403; 
Pellerin et al. 2002). Therefore, we consider it unlikely 
that Star 3 contributes more than 10% to the supernova 
continuum at A > 1350 A. 

We propose that this continuum is predominantly H 
I 2-photon (2s 2 Si/2 - Is 2 Si/ 2 ) emission. We fit the 

1350 < A < 1720 A spectrum with the analytic formula 



of iNussbaumer fc Schmutzl (|1984D . 

■ h = ^ A ^ N ^s (1) 

where J\ is the emissivity and N 2s 2 $ is the H I column 
density in the 2s 2 «S' state. This function rises sharply 
from 1216 A to a maximum at A ~ 1420 A, and then 
declines slowly to the red. The 2-photon spectral model 
is shown in orange in Figure 6. Based on our near-UV 
images, we estimate that the circumstellar ring fills ap- 
proximately 30 ± 20% of the inner 2" of the COS aper- 
ture where the transmission is > 0.6. Assuming this fill- 
ing fraction, the continuum level determines the column 
density, and we find logi (V 2s 2 S (H)) = 21.34^°;^ cm" 2 . 

The total integrated (1216 A - 6 /im) H I 2-photon flux 
is 9.1 ± 0.6 x 10~ x photons cm~ 2 s _1 , and the observed 
ratio of Lya line emission to 2-photon is F(Lya)/F(2s) 
= 1.96 ± 0.23. Since this ratio is consistent with the 
2p/2s ratio of 2.1 (observed as F(Lya)/F (2s)) expected 
for recombination at T » 10 4 K (|Spitzerlll978h and the 
(Is 2p)/(ls -> 2s) ratio of 2.05 -2.10 exp ected for 
excitation by thermal electrons (| C allawavll 1 9 8 81) , we infer 
that it is likely that the Lya and the 2-photon emission 
come from the same source. 

An important constraint on the origin of the far-UV 
continuum is the density of the emitting region. H I 
2-photon emission can be suppressed by a factor [1 + 
Up/ucr]" 1 due to collisional depopulation of the 2s state, 
where the critical density is given by n cr — A2~ t /C sp and 
n p is the density of protons. The (2s —¥ 2p) collision rate 
coefficient is C sp — 5.31 x 10~ 4 cm 3 s _1 (Seaton 1955; 
see also Dennison et al. 2005 for a discussion of 2s and 
2p level populations in H II regions) an d A 2l ~ 8.23 s _1 
for the 2-photon decay (Klarsfeld 1969), therefore n cr = 
1.5 x 10 4 cm^ 3 . 

Based on the observed 2-photon flux, we can predict 
the total amount of associated Ha if this emission is 
predominantly recombination. In Case B recombina- 
tion, the rate of populatio n of the 2s state is 3.3 time s 
the rate of emission of Ha (|Osterbrock &: Ferlandl l2006) . 
The expected Ha in this scenario is « 2.8 x 10~ 4 pho- 
tons cm~ 2 s . We have made an estimate of the total 
Ha flux from the circumstellar ring by analyzing STIS 
G750M spectra obtained in 2009 October (HST expo- 
sure IDs OB7I200A0 - OB7I200D0), making a correc- 
tion for differential reddening between Ha and the 2- 
photon emission. The total Ha emission from the ring is 
w 2.5 x 10 _1 photons cm -2 s _1 , very close to the value 
expected if the 2-photo n emission was produced b y re- 
combination. However, iGroningsson etaTl (|2008bn find 
that the Ha flux from the ring is dominated by high- 
density (> 4 x 10 6 cm~ 3 ) shocked material around 5700 
days after the explosion (October 2002). 2-photon emis- 
sion could be highly suppressed from such gas, which 
has n p /n cr > 300. 2-photon emission from lower den- 
sity (< 5 x 10 3 cm" 3 ; Matilla et al. 2010) gas in the 
unshocked ring would not be suppressed, but such gas 
can only account for ~ 20% of the observed 2-photon 
emission. 

Because we observe no comparable source of Ha emis- 
sion, we conclude that the majority of the 2-photon emis- 
sion is produced not by recombination, but by thermal 
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(T » 10 4 K) electron impact excitation of the 2s 2 S 
state of neutral hydrogen atoms in the outer ejecta, 
near the reverse shock front. This low-density gas 
{n H ~ 100 cm 3 ; Smith et al. 2005; Heng et al. 2006) is 
heated by X-rays emitted from the shocked gas near the 
hotspots. In particular, the energy deposition of the soft 
X-ray/EUV photons from the shocked hotspots will be 
concentrated in a layer near the reverse shock (Fransson 
et al. 2011). When the ionization fraction is > 3 x 10~ 2 , 
the majority of the X-ray energy heats t he gas through 
Coulomb stopping of fast photoelectrons (|Xu et al.lll992t 
iKozma fc Fransson|[i992l ). The primary coolants for this 
gas are Lya and 2-photon emission. 

Is it reasonable to expect this level of X-ray heat- 
ing of gas near the reverse shock? We calculate the 
total X-ray flux from the circum stellar ring us i ng th e 
two-component model spectrum of iZhekov et al.l (|2006[) . 
scaled to the total 0.5 - 2.0 keV luminosity observed by 
Chandra near day 8000 (£(0.5 - 2.0 keV) « 1.5 x 10 36 
erg s _1 ; Racusin et al. 2009). The total 2-photon lu- 
minosity is L(2s) ~ 2.2 x 10 36 erg s -1 , therefore the 
0.5 - 2.0 keV X-ray flux is insufficient to power the far- 
UV continuum. However, the total shock luminosity is 
most likely dominated by emission in the soft X-ray/EUV 
band (0.01 - 0.5 keV) that is attenuated by the neutral 
hydrogen in the interstellar media of the Milky Way and 
LMC (Fransson et al. 2011). The luminosity in this band 
inferred from the model by Zhekov et al. is L(0.01 - 0.5 
keV) « 3 x 10 38 erg s -1 . Assuming that roughly half of 
this emission intersects the outer ejecta, we find that ap- 
proximately 1.4% of the soft X-ray/EUV luminosity from 
the shocked ring must be reprocessed into H I 2-photon 
emission. 

4.1.2. H I Line Emission 

Figures 4 and 5 show the broad emission from H I 
Lya. This emission was first observed in HST- 
STIS spectra 10.25 years a fter the SN 1987A explo- 
sion ([SonncbomeOD [1998) . Subsequent studies have 
explor ed the reverse shock Lya a nd Ha emission in more 
detail (Michael et al. 1998, 2003), including its brighten- 
ing and deceleration (Heng et al. 2006; see also Smith et 
al. 2005 and Fransson et al. 2011 for discussion of the 
Ha evolution). The combination of increased sensitivity 
and spectral resolution and low instrumental background 
of COS enable us to produce the highest-quality ultravio- 
let velocity profiles of the reverse shock emission to date. 
The velocity distribution of the neutral hydrogen emis- 
sion extends from -13000 - +8000 km s" 1 (Figure 5). 
The velocity maxima are much smaller than the initial 
observations by Sonneborn et al. (1998; ± 20000 km s -1 ) 
and are consistent with the decrease in m aximum pro- 
jected velocity observed from 2004 to 2010 (jFrance et al.l 
2010). Using the limited angular resolution of COS, we 
confirm that emission at negative velocities is concen- 
trated towards the northern side of the circumstellar ring 
while positive velocities are concentrated on the southern 
side of the ring. The total integrated (-12000 - +8000 
km s _1 ) Lya flux inferred from Figure 5 is 1.78 ± 0.18 
photons cm -2 s _1 . 

Previous work noted that the Lya/Ha ratios from 
the reverse shock exceed the 5:1 photon pro duction ra- 
tio expected for a Balmer-dominated shock ([Heng et al.l 




1100 1200 1300 1400 1 500 o 1600 1700 1800 
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Fig. 6.— The combined WNW and ENE spectrum of SN 1987A 
(black) and the spectrum averaged over 0.8 A bins chosen to avoid 
narrow features from hotspot emission and interstellar absorption. 
These data have been corrected for interstellar neutral hydrogen 
absorption (logio(N(H)) = 21.43) and interstellar dust extinction, 
assuming E(B — V) = 0.19 and Ry = 3.1. The broad spectral fea- 
tures are produced by a combination of hydrogen 2-photon emission 
(2s — > Is) and ionic emission from the reverse shock. A theoret- 
ical hydrogen 2-photon spectrum is shown overplotted in orange. 
Reverse shock emission features are identified in green. 

[20061: IHeng fc McCravl[200l . iFrance et al.l (f2010h found 

Lya/Ha ratios > 30 from -8000 2500 km s" 1 and 

> 20 from +3000 - +7000 km s _1 for isolated cuts across 
the northern and southern sides of the reverse shock, re- 
spectively. They attributed the enhanced Lya emission 
to a second source of Lya photons. They argued that 
Lya photons from the hotspots are resonantly scattered 
by onrushing hydrogen with a distribution of velocities 
spanning a width Avh ~ 3000 - 9000 km s^ 1 (unlike 
Lya, Ha is not a resonance line and therefore hotspot 
Ha photons pass through the debris freely). 

While the Lya enhancement at the largest negative 
velocities can likely be explained by this mechanism, 
our preceding discussion of the 2-photon continuum sug- 
gests that the majority of the lower velocity (± 7000 
km s _1 ) Lya excess is attributable to the X-ray heat- 
ing of the outer supernova debris. As described above, 
Lya is the primary coolant for mostly neutral hydrogen- 
rich gas excited by soft X-ray/EUV photoelectrons, and 
the agreement between the observed and theoretical 
F(Lya)/ F(2s) ratio therefore argues that thermal elec- 
tron collisions dominate the production of Lya photons 
in the outer ejecta near the reverse shock front. Resonant 
scattering will cause the newly produced Lya photons to 
emerge preferentially in the outward direction, substan- 
tially favoring the blue-shifted wing as observed. 

4.2. He II X1640 Emission 

The broad emission feature that we attribute to He II 
A1640 from the reverse shock is easily visible above the 
1600 - 1700 A continuum in Figure 6. In Figure 7, we 
compare the He II velocity profile with the red side of the 
N V A1240 (§4.3.1) profile. As neutral or singly ionized 
helium cr osses the reverse sho ck, collisions with thermal 
electrons ([Laming et al.|[l996l ) cause the helium to emit 
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TABLE 2 

Continuum Subtracted SN 1987A Reverse Shock Emission 



Species 







Line Flux 


Flux Uncertainty 




/ i \ 

( A ) 


(km s *■) 


(photons cm * s L ) 


(photons cm ' s L ) 


H I b 


1169.6 - 1208.0 


-12000 - -1500 


1.1 X 10° 


5.9 x 10" 3 


H I 


1222.8 - 1247.1 


1500 - 8000 


4.5 X 10" 1 


4.1 x 10" 3 


N V c 


1197.9 - 1233.9 


-10000 - -1000 


5.7 x 10" 2 


6.2 x 10" 3 


N V 


1247.1 - 1285.5 


2000 - 12000 


8.0 x 10" 2 


8.7 x 10" 4 


N IV] 


1428.3 - 1478.8 


-12000 - -1500 


1.8 x 10~ 3 


4.9 x KT 4 


N IV] d 


1497.0 - 1523.8 


1500 - 8000 






c iv 


1491.9 - 1539.4 


-12000 - -1500 


4.6 x 10~ 3 


5.4 x HT 4 


C IV 


1558.9 - 1584.6 


1500 - 8000 


2.1 x 10~ 3 


4.8 x 10" 4 


He II 


1590.4 - 1633.7 


-10000 - -1000 


6.8 x 10" 3 


7.8 x 10" 4 


He II 


1658.0 - 1705.0 


2000 - 12000 


9.6 x 10~ 3 


8.1 x 10~ 4 



a v = assumed to be the rest wavelength of the stronger lines of the N V and C IV doublets. 

b Corrected for interstellar neutral hydrogen absorption of logio(N(H)) = 21.43, however the inner ± 6 A (± 1500 km s — 1 ) cannot be 
reconstructed. All reverse shock fluxes have been corrected for interstellar dust extinction, assuming E(B — V) = 0.19 and Ry = 3.1. 
c Blue N V flux is overwhelmed by reverse shock Lya emission and interstellar H I absorption. We assume the N V Blue/Red ratio is the 
same as the He II reverse shock emission ( PS 0.7). 

d Red N IV] flux is obscured by the stronger blue emission from C IV. 



the analog to Ha. The velocity distribution of the ion- 
ized helium is expected to be quite different from that of 
neutral hydrogen because the He + ions that have crossed 
the reverse shock will be deflected by turb ulent magnetic 
fields in the shock isotropization zone (Michael ct al. 
f998). The broad He II displays a somewhat asymmetric 
line shape, spanning approximately —9000 — hUOOO km 
s- 1 . 

Interpolating the integrated red and blue reverse shock 
He II A1640 photon fluxes given in Table 2 through 
the narrow emission lineifl, we measure a total flux of 
1.9 x 10 -2 photons cm -2 s _1 . Because this line is the hy- 
drogenic analog to Ha, the F (1640) /F (Ha) ratio should 
give a direct measure of the relative abundance of he- 
lium (by number), assuming full hydrogen and helium 
ionization and that our line identification and contin- 
uum subtraction are correct. We compare our He II 
measurement with the Ha flux from the day ~ 8000 
VLT-UVES observations presented by Fransson et al. 
2011, F obs (Ha) w 1.4 x 10~ 13 erg cm" 2 s _1 . We con- 
vert this value into a total reverse shock Ha by applying 
a factor of 1.86 slit c orrection, a factor of 1.45 redden- 
ing correction at Ha (|Groningsson et ail l2008al) . and a 
factor of 24/22 to account for the relative time of ob- 
servation after the explosion. The integrated Ha photon 
flux in 2011 is w 0.14 photons cm~ 2 s" 1 . The He/H 
abundance ratio in the reverse shock is 0.14 ± 0.06, 
in agreement wi t h the He/H abundance ratio derived 
by iMattila et al.l ([2010) for the circumstellar ring. The 
fact that the He II A1640/Ha ratio reproduces the cir- 
cumstellar ring He/H abundance ratio sugg ests that the 
neutra l hydrogen pre-ionization predicted bv lSmith et al.l 
(120051) is negligible 22 ye ars after the explosion. 

Bor kowski et al.l (|1997[ ) predicted the time evolution of 
the observable (accounting for attenuation due to inter- 
stellar reddening) reverse shock flux in several abundant 
ions, including He II. Their predictions only extend to 
2007, but we can extrapolate their curves to 24 years 
after the explosion for comparison with the COS data. 

3 The total emission from the reverse shock is interpolated 
through the central region lost to narrow-line emission. In the case 
of He II, F§S n = (1 + (3000/22000)) x (F^fed + p RSJlu^ 




Velocity (km s ! ) 



Fig. 7. — 2-photon continuum-subtracted N V and He II in the 
combined WNW + ESE spectra. The blue side of the N V dis- 
tribution is lost under the RS Lya emission. The He II velocity 
profile appears qualitatively similar to that of N V. 

We estimate their prediction for the observed He II flux 
in 2011 to be ~ 1 - 3 x 10~ 3 ph otons cm -2 s _1 from 
Figure 4 of lBorkowski et aT] (|1997f) . Correcting these val- 
ues by a factor of 4 to account for interstellar reddening, 
their prediction for the He II flux is 0.4 - 1.2 x 10~ 2 
photons cm~ 2 s _1 , only a factor of ~ 2 below the He II 
flux observed by COS. Differences between the predicted 
and observed He II fluxes are likely due to assumed elec- 
tron temperatures that are lower than suggested by the 
observations (see the next subsection). 

4.3. Carbon and Nitrogen Emission Lines 

4.3.1. NVX1240 

Prior to the first ST1S observations of SN 1987A, 
Borkowski et al. (1997) predicted that strong reverse 
shock emission from Li- like N V A 1240 would be de- 
tectable. This emission was not apparent, however, 
in the first deep far-UV STIS spectra presen ted by 
iSonneborn et al.l ([19981 ) and iMichael et al.l (|1998l) . In re- 
cent (January 2010) STIS observations, we tentatively 
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detected broad, redshifted N V emission (Figure 3 of 
France et al. 2010), but low S/N precluded a detailed 
analysis. Now, in our COS observations, we unambigu- 
ously detect this emission. While the blue wing of the 
line is lost under the bright Lya emission, we observe the 
complete red wing of the N V A1240 velocity profile. In 
Figure 7, we compare the line profile of N V with that of 
He II, observing that the red wings of the two profiles are 
qualitatively similar. Both ions present a boxy line pro- 
file, with a flat top and fall off between +9000 - +10000 
km s^ 1 . The N V profile extends to at least +14000 km 
s _1 , where the data are contaminated by the geocoronal 
O I triplet. 

At the relevant energy scales of nonradiative super- 
nova shocks, v s > 10 3 km s _1 , the excitation cross sec- 
tions for Li-like species (e.g., C 3+ , N 4+ , 5+ ) by col- 
lisions with protons and heavier ions are considerably 
larger than the cross-section fo r ionization to their He - 
like stage (jLaming et all 119961: iBorkowski et al.l 119971 ) . 
Therefore, these species may emit many (several hun- 
dred) line photons (predominantly C IV A 1550, N V A 
1240, O VI A 1032 for the ions listed above) for every 
atom that crosses the reverse shock front before becom- 
ing ionized. The total ratio of N V A1240 photons emit- 
ted to Ha is given by 



cm" 3 dSmith et all 120051: iHeng et al.1 120061). and n He 



F(NV) 



1 x N 



#(1240) 



F(Ha) 0.2 x H R{N A + N 5 +) 



(2) 



where 0.2 is the number of Ha photons e mitted per neu- 
tral hydrogen atom crossing the shock (Michael et al. 
|2003f ), xn/xh is the nitrogen abundance ratio, i?(1240) 
is the rate of electron and ion impacts that result in N V 
emission, and R(N 4+ — > N 5+ ) is the rate of ionization 
to N 5+ . The relative nitrogen abun dance in the circum - 
stellar ring is 2.8 (± 1.1) x 10" 4 (jMattila et al.ll2010f ). 
Given the similarity of the red side of the He II and N V 
profiles, we compute the total N V reverse shock pho- 
ton flux by assuming that the relative red/blue contri- 
bution is the same as He II (Table 2) and interpolating 
through the narrow emission lines, as described above. 
Under these assumptions, the observed F(N V)/F(Ha) 
ratio is » 1.15. Folding in the nitrogen abundance, 
we can therefore make a rough estimate of the ratio 
i?(1240)/i?(iV 4+ -)■ N 5+ ) ~ 850±||g required to explain 
the observed line ratio. 

i?(1240) can be evaluated analytically, 



#(1240) 



(cr e v s ) + n H (a p v s ) + n He (a a v s ) (3) 



where the cross-sections, <J e ,p,a, are for excitation of N V 
A1240 by electrons, protons, and alpha particles, respec- 
tively. We simplify the averages over velocity distribution 
by evaluating the expression at a single shock velocity, v s 
= 10 4 km s . Cross-sections for the pro ton and alpha 
partic le collisions are taken from Table 2 of lLaming et al.1 
(1996) for proton energies (E p — 1/2 m p v 2 ) of 522 keV 
(we assume the 544 keV values) and alpha particle en- 
ergies (E a = 1/2 (Am p )v 2 s ) of 2088 keV (we assume 1.2 
times the 2720 keV values). The electron impact ex- 
citation cross-sections are not as well determined, but 
we expect them to be similar to those of the protons 
at high shock velocity. Therefore, we estimate a e ,p,u 
= [4.31, 4.31, 14.4] x IP" 17 cm" 2 , respec tively. Tak- 
ing n e = 450 cm -3 (jBorkowski et al.1 119971 ). n H w 100 



0.17 nn (jMattila et al.ll2010D . we compute a total N V ex- 
citation rate for v s = 10 4 km s" 1 ; i?(1240) = 2.6 x 10~ 5 
s" 1 . 

Ionization rates by baryons are expected to be small 
compared to those by electrons due to th e high energy 
of po st-shock protons and heavier species (|Laming et al.1 
1996), therefore we neglect their contribution to the 
N 4+ ionization rates. Cross-sections for the direct 
ionization to N 5+ by thermal electrons can be calcu- 
lat ed as a function of electron energy from Equation 1 
of lArnaud fc Rothenfluel ([1985). Using coefficients for 
the Li-sequence (their Table 1), we calculate an electron 
ionization cross-section of a, 



NV 



at 



107 eV0. The direct ionization rate, R(N 



4+ 



) = 1.5 x 10~ 7 s" 1 . This gives 



E e = 
N 5 +) 

R(1240)/R(N 4+ -> N 5+ ) ~ 170, which is outside the 
range allowed by the observations. 

If, on the other hand, there is partial equili- 
bratio n between the electron and proton distribu- 
tions (jCargill fc Papadopouloslll988HLaming et al.lll996l : 
[He^[2010h . then T' e 



PeqT p , where 

3 iim p v1 
16 



(4) 



and (3 eq is the equilibration factor, which has a maximum 
value of 1. fj, is the mean particle weight, fi ~ 0.55 (for 
Xh = 0.85, xhc — 0.14, and xz = 0.01). In this case E e 
= 1/2 m e (2/csT e / m e ) — kBT e . Using this energy to 
calculate the ionization cross-sections, we find that for 
f3 eq = 0.14 - 0.35 (T' e = 1.7 - 4.4 x 10 s K), the elec- 
tron ionization cross-sections cr^X ~ 9-7 - 4.2 x 10~ 20 
cm 2 , and the ionization rates of R(N 4+ -> iV 5+ ) = 4.3 - 
1.9 x 10~ 8 s _1 . This gives R(1240) / R(N 4+ -> iV 5 +) 
~ 610 - 1380, approximately the range required to ac- 
count for the observed F(N V)/F(Ha) rati© v s = 10 4 
km s _1 is used as the fiducial velocity in our calculations, 
but the general conclusions hold for a range of possible 
shock velocities. Computing i?(1240)/i?(7V 4 + -> iV 5+ ) 
for velocities (5 - 12) x 10 3 km s _1 , we find that /3 eg 
decreases by ~ 50% across this range, but (3 eq > 0.1 is 
required to explain the data for all velocities considered. 

It seems clear that partial ion-electron equilibration 
is required to explain the observation of strong reverse 
shock N V emission from SN 1987A. However, values 
of B eo = 0.14 - 0.35 are greater than those favored 
by ILaming et al.l (| 19961 ) from their non-radiative shock 
modeling of SN1006. Moreover, equilibration values of 
Peq > 0.1 for the high shock velocity (v s ~ 10 4 km 
s^ 1 ) in SN 1987A are inconsistent with some recent re- 
sults on electron-ion equilibration in collisionless shocks. 
Ghavamian et al. (2007) discuss a sample of older (~ 10 3 
yr) supernova remnants whose electron-to-proton tem- 
perature ratio is inversely proportional to the square of 



the shock velocity, j3, 



eq 



however, this relation 



5 We note that if we assume that the ionization cross-sections 
scale as cr^^ oc T e ~ 0,5 in the high-T e limit, then the ion-electron 
equilibration factor could be a factor of 2 - 3 higher, approaching 
the limit of complete equilibration. 
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has n ot been confirm ed in subsequent wor k ( Hcldc r~et al.l 
1201 If) . Additionally. Ivan Adelsberg et al.l (|2008T ) present 
evidence suggesting an increasing (5 eq for v s > 2000 km 
s . Finally, we note that very high electron temper- 
atures (T e ~ 10 9 K) were inferred for the high veloc- 
ity ejecta in SN 1993J (|Fransson et alJll996f) . suggesting 
that electron heating may be efficient in high velocity 
environments (see §5.3). 

We also compare the obs e rved N V flux with that pre- 
dicted by [Borkows ki et al.l (jl997T ) . Extrapolating their 
prediction to the time of our observations, we would have 
expected (1.5 - 2.0) x 10~ 13 erg cm -2 s _1 in 2011, or 
0.06 - 0.08 photons cm~ 2 s _1 after including a factor of 
6.3 to account for interstellar reddening. We see that, 
similar to the prediction for He II, these values are a fac- 
tor of ~ 2 below their observed 2011 values (F(N V) = 
0.16 photons cm -2 s _1 ). Overall, we consider the predic- 
tion made 14 years prior to the observation to be remark- 
ably good. The slight differences between the predicted 
and observed N V flux can most likely be attributed to 
their choice of a value for /3 eq that is lower than we de- 
termine from the data. Lower electron energies increase 
the ionization cross-section and hence R(N 4+ — > N 5+ ), 
thereby reducing the total number of N V photons emit- 
ted. 



4.3.2. NIV]X1486 

Figure 6 shows that there is considerable substructure 
on the continuum in addition to the strong reverse shock 
emission from HI, He II, and N V. The weakest of these 
features spans ~ 1430 - 1500 A. We attribute this emis- 
sion to N IV] A1486 from the reverse shock. The N IV] 
emission is highly asymmetric about the rest velocity. 
The blue-shifted component extends to roughly —8000 
km s _1 , while the red side of the profile is at the noise 
level by +3000 km s _1 . Additionally, the blue wing of 
the reverse shock C IV profile, discussed in the next sec- 
tion, overwhelms any additional red flux. Again making 
the assumption that the red/blue N V ratio is the same 
as that for He II, and assuming that the red side of the 
N IV] profile contributes an additional 20% to the total 
flux, we find that the total F(N V)/F(N IV]) ratio is 
~ 72 ± 17. 



4.3.3. CIVX1550 

Figure 6 also shows C IV A1550 from the reverse shock 
covering ~ 1500 - 1580 A, with emission from —12000 - 
— 1500 km s" 1 having 2.2 times the total integrated pho- 
ton flux as the +1500 - +8000 km s" 1 interval (Table 
2). The red side of the velocity profile is only marginally 
above the noise level at vciv > +3000 km s _1 . We 
consider the relative fluxes of C IV and N V here, and 
in the next section we will discuss the relative velocity 
distributions of the reverse shock species. 

Interpolating the spectrum through the narrow emis- 
sion lines, we measure an integrated broad line photon 
flux (-12000 - +8000 km s" 1 ) of F(C IV) = 7.7 x 10" 3 
photons cm -2 s _1 . With the previously-noted assump- 
tions about the N V flux in mind, we find a flux ratio of 
F(N V)/F(C IV) » 20 ± 3. The ratio of N V A1240/C IV 
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Fig. 8. — 2-photon continuum subtracted N IV], C IV, and He II 
in the combined WNW + ESE spectra. He II and C IV have 
been offset by 4.0 and 2.0 photons cm -2 s — 1 A -1 , respectively, for 
display purposes. The zero-flux levels for all three ions are shown 
as the dashed gray line. 



A1550 is given by an expression similar to equation (2): 



F(NV) 
F(CIV) 



XN 

xc 



#(1240) 



R(N 4 + -> N 5 +) 



#(1550) 



R(C 3 + -> C 4 +) 
(5) 

The ratio of the quantities in square brackets is very close 
to 0.9 for v s = (5 - 12) x 10 3 km s _1 (the ratio is near 
unity because both are Li- like ions). That implies that 
the abundance ratio of nitrogen/carbon atoms crossing 
the reverse shock is xn/xc ~ 22 ± 3, greater than the 
ratio xn/xc ~ 8.5 ± 3.5 inferred from spectroscopic ob- 
servations of the equatorial ring (Lundqvist & Fransson 
1996; Mattila et al. 2010). 

We note that while our estimate of the flux ratio 
N V A1240/C IV A1550 is somewhat uncertain because 
our choice of the blue side of the N V profile is specula- 
tive, the integrated flux in the red side of the N V profile 
alone is more than 10 times greater than the total C IV 
flux. Therefore, we are confident that the N V/C IV 
ratio is enhanced well beyond what can be attributed to 
the circumstellar ring abundances of the two species. We 
return to this point in Section 5.4. 

5. DISCUSSION 
5.1. Relative Velocity Distributions 

In Sections 4.2 and 4.3, we describe the velocity distri- 
butions of He II, N V, N IV], and C IV. The red side of 
the N V velocity profile is compared with He II in Fig- 
ure 7 and a comparison of the He II, N IV], and C IV 
velocity profiles is shown in Figure 8. We observe quali- 
tatively different velocity profiles from He II, N IV], and 
C IV. He II displays a mildly asymmetric profile span- 
ning 9000 - +11000 km s _1 , while N IV] and C IV are 

strongly weighted towards the blue, spanning ~ —(9000 - 
8000) — h3000 km s _1 . Peculiar geometric projection of 
the velocity can be ruled out because all three ions were 
observed with the G160M instrument mode, with a con- 
stant spacecraft roll orientation. 

In nonradiative shocks, He II is thought to be excited 
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by collisions with electrons, while the Li-li ke species (N V 
and C IV) are excited by ion-collisions (|Laming et al.l 
1996). As we have shown above, the electron and ion dis- 
tributions have likely not equilibrated in the SN 1987A 
debris, therefore it may be that the relative distributions 
of the impacting species is reflected in the velocity pro- 
files of the excited ions. One might also speculate that 
if ion collisions dominate the emission from metals, the 
preferentially blue N IV] and C IV profiles would reflect 
an asymmetric spatial distribution of the protons and 
ions. The N IV] and C IV profiles also show spatial vari- 
ations between the WNW and ENE spectra. In contrast, 
the He II profile, which is primarily excited by electrons, 
is more spatially and spectrally symmetric. This may 
indicate a more uniform distribution of electrons. With 
these possibilities in mind, one is tempted to separate 
the metals from the helium, however we observe N V ex- 
tending to red velocities as high as +10000 km s _1 before 
declining. Modeling of both the fluxes and velocity dis- 
tributions of ions in the SN 1987A reverse shock would 
be very useful for our understanding of the observed ve- 
locity profiles. 

5.2. Excess Emission at A < 1350 A 

The high throughput and low instrumental background 
of COS have allowed us to study the ultraviolet con- 
tinuum of SN 1987A for the first time. Similarly, we 
have presented the first conclusive observation of ele- 
ments heavier than hydrogen in the reverse shock. Fig- 
ures 3 and 6 show that there is additional ultraviolet 
continuum at A < 1350 A, above which we propose that 
hydrogen 2-photon emission can adequately explain the 
data. As we discuss in §4.1.2, the spectral contamina- 
tion from Star 3 is inversely proportional to wavelength, 
as the cross-dispersion astigmatism height is largest at 
the shortest wavelengths. This could account for some, 
but probably not all, of the bluest continuum emission. 
Instead, we propose four speculative possibilities for this 
emission. He + has a 2p — > Is Lya line at 304 A, and 
will emit a 2s — > Is 2-photon continuum which peaks at 
A « 355 A, and declines to the red. We estimate that this 
emission contributes « 10% of the 1150 A continuum. 

There may be an additional Lya component at 
very high velocities, although this would require 
VLya > 30 x 10 3 km s" 1 (> 0.1 c) to account for 
the reddest observed emission. Another possibility is 
that low-ionization line emission from the core, as op- 
posed to higher-ionization emission from atoms cross- 
ing the re verse shock, contr i butes to the A < 1350A 
emission. IJerkstrand et al.l ()2011[ ) describe 2-photon 
emission, scattered by low-ionization, high-opacity metal 
lines in the core, contributing to the observed spectrum 
at ~ 8 years following the explosion. The observed ex- 
cess begins at roughly the wavelength of the O I] A1356 A 
emission line. If an appreciable amount of low-ionization 
metal emission escapes from the near side of the core , it 
may be responsible for a portion of the observed flux. 

Perhaps the most likely scenario is one in which sev- 
eral minor reverse shock species contribute to the ob- 
served continuum between 1140 - 1350 A. We have seen 
that reverse shock emission from Be-like ions (e.g. NIV] 
A1486) are present. Therefore, emission from lines of 



C III A1175 and O V A1371 may also contribute^. 

5.3. Partial Ion-Electron Equilibration in SN 1981 A: 
Evidence for the Cross- Shock Potential? 

In Section 4.3.1, we found that large ion-electron tem- 
perature equilibration ratios (f3 eq » 0.14 - 0.35) are re- 
quired to explain the high N V A 1240 flux observed 
in the spectrum of SN 1987A. IGhavamian et al.l (j2007[ ) 
predict equilibration ratios several times less than this 
for shock velocities observed in SN 1987A (v s ~ 10 4 
km s _1 ). Therefore, we consider additional electron 
heating mechanisms that could explain the observed 
N V emission. One mechanism that has been proposed 
for collisionless shocks (usually assumed to be quasi- 
perpendicular) at low Alfve n mach numbers (M a) is 
the cross-shock potential fe.g lBalikhin et aLlll993l) . But 
the cross-shock potential is expected to become less im- 
portant as Ma increases because the plasma resistivity 
causes i nsufficient dissipation and the shock becomes un- 
steady ()Edmiston fc Kennell 119841 ) . Electrons undergo- 
ing E x B drift across the shock front are not guaran- 
teed to see the same cross-shock potential. They may 
lose energy as well as gain it, and this limits the electron 
heating. 

However, it has also been suggested that at Alfven 
Mach numbers appropriate to shocks in supernova rem- 
nants, the shock transition becomes thin. The shock 
transition can develop a length scale less than the 
electron convective gyroradius, thus eliminating any 
E x B drift within the shock. This onsets at Ma 
= ^Jm^Jm^ and electron heatin g by the cross shock 
potential aga i n bec omes viable (jBalikhin et al.l 119931 : 
iGedalin etaLl 120081 ). Typical SNR shocks have M A 
much higher than this, neglecting any modification of 
the preshock medium. Yet the predicted heating is not 
observed. Most likely, Ma at these shocks is much lower 
than expected due to amplification of the ma gnetic field 
by co smic rays streaming ahead of the shock ( BelH l2004l 
2005). Electron heating in these systems occurs via 
plasma waves excited upstream by shock reflected ions 
or cosmic rays. At higher velocity shocks, or for those 
likely to have weak upstream magnetic field due to their 
environment and/or insignificant cosmic ray populations, 

Ma > \J rm/rrie and the cross-shock potential may be ca- 
pable of heating electrons to the required energy There- 
fore the energetic environment of SN 1987A, especially at 
its reverse shock, may enable additional electron heating 
mechanisms that do not contribute in lower Ma systems. 

5.4. The N V/C IV Ratio 

The final unresolved issue is the large N V/C IV ra- 
tio in the reverse shock emission. Carbon and nitro- 
gen abundance ratios in t he circumstellar ring suggest 
this number should be w 8 ([Lundqvist fc Franssoij|l996t 
iMattila et al.l I2010D , however the observed flux ratios 
indicate N/C ss 22. It may be that the N/C abun- 
dance ratio in the outer envelope of the progenitor was 
stratified prior to the ejection of the circumstellar rings, 
and that we are now seeing the first observational ev- 
idence of that stratification. A second possibility to 

6 O V is the next ionization stage down from the Li-like O VI 
A1032 line which we do not detect, but is presumably in the reverse 
shock spectrum at shorter wavelengths. 



COS Observations of SN 1987A 



11 



account for this discrepancy is that ongoing thermonu- 
clear processing continued to convert C to N in the su- 
pernova progenitor following the circumstellar ring ejec- 
tion. The CNO bi-cycle will enrich the 14 N abundance 
at the expense of the abundances of 12 C and 16 and, 
in equilibrium, will c onvert almost all of the primor- 
dial C and O into N (ICauehlan fc Fowled fl962h . CNO 
processing has been invoked to explain the fact that 
the observed nitrogen abundance in the circumstellar 
ring is elevated by a factor ~ 10 over its value in the 
LMC (jFransson et al.l fl989) . Our observation that the 
He abundance does not change between the ring and 
reverse shock is qu alitatively consistent with th e 14E1 
model presented by iShigevama fc Nomotol (|1990() . They 
show that the high-velocity material we observe crossing 
the reverse shock front is only a small fraction of the to- 
tal ejected mass, and does not probe deep enough into 
the interior of the ejecta to observe significantly elevated 
He abundances. 

If CNO processing continued near the stellar surface 
following the ejection of the circumstellar rings, it could 
have in principle converted most of the remaining C 
and O abundanc es seen in the equatorial ring into N. 
iHeng et all (|2008l ) note reduced oxygen abundances, pos- 
sibly related to the high N V/C IV ratio observed in 
our observations. For this explanation to be viable, the 
timescale to reach equilibrium in the CNO cycle must 
be < 20,000 years, the interval since the ejection of the 
equatorial ring. This condition will be met if the tem- 
perature of the shell where CNO burning takes place is 
T > 3.5 x 10 7 K (Caughlan and Fowler 1962, Table 4). 
This effect could not only account for the high N/C ra- 
tio in the gas crossing the reverse shock, it could also 
increase the N/H ratio by an additional factor of ~ 2. 
This would lower the required /3 eq derived in §4.3.1. 

6. CONCLUSIONS 

We have presented an analysis of deep HST-COS spec- 
troscopy of SN 1987A. Below, we summarize our primary 
results. 

1. We observe narrow lines from shocked gas in the 
circumstellar ring, broad emission lines from the 
reverse shock, and a strong detection of ultraviolet 
continuum emission. Several of the emission line 
species and the far-UV continuum are conclusively 
detected for the first time. 

2. The asymmetry in Lya profile and its enhancement 
relative to Ha suggests that most of the Lya emis- 
sion is a result of the illumination of the outer su- 
pernova debris by soft {E < 0.5 keV) X-rays emit- 
ted by the shocked circumstellar ring. 

3. The A > 1350 A continuum may be described 
by hydrogen 2-photon emission originating in the 
outer ejecta, near the reverse shock front. 

4. We present resolved velocity profiles of He II, C IV, 
N IV], and N V from the reverse shock. A com- 
parison of the velocity integrated He II and Ha 
velocity profile indicates a He/H abundance ratio 
of 0.14 ± 0.06. In order to reproduce the observed 



F(N V) /F(Ka) line ratio, partial ion-electron equi- 
libration is required, T e = (3 eq T p , with j3 eq « 0.14 - 
0.35. Large values of f3 eq may be explained by addi- 
tional electron heating by the cross-shock potential. 

5. The velocity profiles of C IV and N IV] are signif- 
icantly different from that of He II, which may be 
related to the different excitation processes (ion vs. 
electron collisions) for the different species. 

6. We observe additional continuum emission at 
A < 1350 A that is not readily explained. We fa- 
vor a scenario where several weaker emission lines 
contribute, but spectral overlap with the much 
stronger Lya and N V emission profiles prevents 
an unambiguous interpretation. 

7. Finally, we observe that the C IV reverse shock 
emission is weaker than would be expected assum- 
ing circumstellar ring abundances. This may be ex- 
plained by chemical stratification in the outer enve- 
lope of the progenitor, and may indicate additional 
CNO processing between the period of circumstel- 
lar ring ejection and the supernova explosion. 
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